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The  tilt  angle  of  a  solar  energy  system  is  one  of  the  important  parameters  for  capturing  maximum  solar 
radiation  falling  on  the  solar  panels.  This  angle  is  site  specific  as  it  depends  on  the  daily,  monthly  and 
yearly  path  of  the  sun.  The  accurate  determination  of  the  optimum  tilt  angle  for  the  location  of  interest  is 
essential  for  maximum  energy  production  by  the  system.  A  number  of  methods  have  been  used  for 
determining  the  tilt  angle  at  different  locations  worldwide.  Keeping  in  view  the  relevance  of  the 
optimum  tilt  angle  in  energy  production  and  reducing  the  cost  of  solar  energy  systems,  the  present  study 
has  been  undertaken.  This  paper  provides  the  update  status  of  research  and  applications  of  various 
methods  for  determining  solar  panel  tilt  angle  using  different  optimization  techniques.  The  study  shows 
that  for  maximum  energy  gain,  the  optimum  tilt  angle  for  solar  systems  must  be  determined  accurately 
for  each  location.  The  review  will  be  useful  for  designers  and  researchers  to  select  suitable  methodology 
for  determining  optimal  tilt  angle  for  solar  systems  at  any  site. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


Contents 


1.  Introduction . 503 

2.  Solar  radiation  modeling  on  sloped  surfaces . 504 

2.1.  Solar  time . 504 

2.2.  Solar  geometry . 504 

2.3.  Global,  beam,  diffuse  irradiance  and  the  sky  radiance  distribution . 504 

2.4.  Solar  radiation  on  sloped  surfaces . 504 

2.4.1.  Isotropic  model . 505 

2.4.2.  Anisotropic  model . 506 

3.  Optimum  tilt  angle  determination . 506 

3.1.  Total  solar  radiation  on  tilted  surface . 506 

3.2.  Tilt  angle-latitude  relations . 507 

3.3.  Tilt  angle  optimization  using  GA,  SA  and  PSO  techniques . 508 

3.4.  Tilt  angle  optimization  using  ANN  techniques . 508 

4.  Overview  of  optimum  tilt  angle  determination  methods . 508 

5.  Conclusions .  511 

References .  511 


1.  Introduction 

The  design  and  development  of  solar  energy  systems  require 
the  knowledge  of  variation  and  maximum  utilization  of  solar 
radiation  falling  on  it.  The  measured  solar  radiation  data  are  not 
available  for  large  number  of  sites,  so  it  has  to  be  estimated.  A 
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number  of  techniques  are  used  for  the  estimation  of  solar  radia¬ 
tion  on  horizontal  surface  [1-3].  The  amount  of  solar  radiation 
received  by  a  solar  photovoltaic  [PV]  panel  or  a  solar  thermal 
collector  is  mainly  affected  by  its  orientation  and  tilt  angle  [4,5], 
The  solar  panels  are  generally  oriented  toward  the  equator,  in  the 
northern  hemisphere  oriented  toward  south  and  in  the  southern 
hemisphere  toward  north.  However,  the  solar  radiation  is  site 
specific  with  diurnal,  monthly,  seasonal  and  yearly  variations;  as 
such  the  optimum  tilt  angle  for  capturing  maximum  solar  radia¬ 
tion  will  also  vary  for  every  location.  This  study  provides  an  update 
status  of  research  and  applications  of  various  computational 
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Nomenclature 

^ext 

extra-terrestrial  daily  radiation  incident  on  a  horizon- 

tal  surface,  Wh/m2/day 

ft  opt 

optimum  tilt  angle 

he 

solar  constant  1367  W/m2 

GA 

genetic  algorithm 

h 

total  radiation  on  tilted  surface 

H 

monthly  average  daily  total  radiation  on  horizontal 

n 

nth  day  of  the  year 

surface 

PSO 

Particle  Swarm  Optimization 

Hb 

daily  beam  radiation  incident  on  an  inclined  surface, 

PV 

photo  voltaic 

Wh/m2/day 

Rb 

ratio  of  average  daily  beam  radiation  incident  on  an 

Hb 

monthly  average  daily  beam  radiation 

inclined  surface  to  that  on  a  horizontal  surface 

Hd 

daily  sky-diffuse  radiation  incident  on  an  inclined 

R 

ratio  of  total  radiation  on  tilted  surface  to  total  radia- 

surface,  Wh/m2/day 

tion  on  horizontal  surface 

Hd 

daily  diffuse  radiation  incident  on  a  horizontal  surface, 

R 

ratio  of  monthly  average  daily  radiation  on  tilted 

Wh/m2/day 

surface  to  monthly  average  daily  total  radiation  on 

Hd 

monthly  average  daily  diffuse  radiation 

horizontal  surface  ( HT/H ) 

"g 

daily  global  radiation  incident  on  a  horizontal  surface, 

Rb 

ratio  of  average  daily  beam  radiation  on  tilted  surface 

Wh/m2/day 

to  horizontal  surface  for  a  month 

Hr 

daily  ground  reflected  radiation  incident  on  an 

Rd 

ratio  of  average  daily  diffuse  radiation  incident  on  an 

inclined  surface,  Wh/m2/day 

inclined  surface  to  that  on  a  horizontal  surface 

Hs 

sky  diffuse  radiation 

Wsr 

sunrise  hour  angle  (°) 

Ht 

daily  global  radiation  on  a  tilted  surface,  Wh/m2/day 

wss 

sunset  hour  angle  (°) 

Ht 

monthly  average  daily  radiation  on  tilted  surface 

p 

surface  slope  from  the  horizontal 

I 

total  irradiation  on  horizontal  surface  (J/m2) 

s 

declination  angle  (°) 

Id 

diffuse  radiation  on  horizontal  surface  (J/m2) 

4> 

latitude  (rj 

ld,T 

diffuse  radiation  on  tilted  surface 

p 

ground  albedo 

Y 

surface  azimuth  angle  (°) 

methods  for  determining  solar  panel  tilt  angle  using  different 
optimization  techniques.  This  paper  is  organized  as  follows: 
Section  2  describes  the  methods  to  determine  solar  radiation  on 
sloped  surfaces.  Section  3  provides  details  of  optimum  tilt  angle 
determination.  An  overview  of  optimum  tilt  angle  determination 
methods  is  given  in  Section  4.  The  conclusions  of  the  study  are 
given  in  Section  5. 

2.  Solar  radiation  modeling  on  sloped  surfaces 

In  this  section,  basic  solar  radiation  parameters  and  models  to 
determine  solar  radiation  on  sloped  surfaces  are  described.  The 
solar  radiation  incident  outside  the  earth's  atmosphere,  called 
extraterrestrial  radiation  for  a  day  of  the  year,  n  is  given  as 


the  surface  from  the  local  meridian,  with  zero  due  south,  east 
negative  and  west  positive,  is  called  surface  azimuth  angle  (p); 
-180°<y<180°.  The  hour  angle  (w)  is  the  angular  displacement  of 
the  sun  east  or  west  of  the  local  meridian  due  to  the  rotation  of  the 
earth  on  its  axis  at  15°  per  hour  with  morning  negative  and 
afternoon  positive. 


2.3.  Global,  beam,  diffuse  irradiance  and  the  sky  radiance 
distribution 

The  distribution  of  solar  radiance  over  the  sky  (Fig.  2)  consists 
of  three  components:  isotropic  dome,  circumsolar  and  horizon 
brightening  [7].  The  circumsolar  and  horizon  brightening  are 
assumed  to  be  concentrated  at  the  center  of  the  sun  and  at  the 


J/JC[1.0+0.333cos(360n/365)]  1 

leM=  \  /jc[1  .00011 0  +  0.03421  lcos((n-l)(360/365))+0.00128sin((n-l)(360/365))+ 0.00071 9cos  2((n-l)(360/365))  +  0.000077sin  2((n-l)(360/365))  J 

2.1.  Solar  time  (1) 


The  solar  time  is  given  as  follows: 

Solar  time  =  standard  time+4(Lst-Ltoc)+£  (2) 

E  =  229.2(0.000075  +  0.001 868cos  B-0.032077 

sin  B-0.014615cos  2B-0.04089sin  2 B)  (3) 

where  B=(n-l)(360/365),  Lst  is  the  standard  meridian  for  the  local 
time  zone,  Uoc  is  the  longitude  of  the  location. 

2.2.  Solar  geometry 

The  position  of  the  sun  in  the  sky  is  given  by  zenith  angle  ( 9Z ) 
(Fig.  1)  [6].  The  angular  position  of  the  sun  at  solar  noon  is  called 
declination  (5);  -23.45  <<5<23.45  .  The  angle  between  the  plane  of 
the  surface  and  horizontal  is  called  slope  (/?);  0°</Kl80°.  The 
deviation  of  the  projection  on  a  horizontal  plane  of  the  normal  to 


horizon.  The  horizon  brightening  for  clear  skies  is  the  highest  at 
the  horizon  and  its  intensity  decreases  away  from  the  horizon  and 
for  overcast  skies.  The  sky  radiance  increases  instead  of  decreasing 
away  from  the  horizon.  The  anisotropic  models  include  diffuse  sky 
radiation  in  the  circumsolar  and  horizontal  brightening  compo¬ 
nents  of  the  solar  radiation.  The  isotropic  models  assume  that  the 
intensity  of  diffuse  sky  radiation  is  uniform  over  the  sky  dome. 

2.4.  Solar  radiation  on  sloped  surfaces 

The  determination  of  solar  radiation  on  sloped  surface  requires 
the  values  of  total  solar  radiation  on  the  horizontal  surface  and  the 
direction  from  which  beam  and  diffuse  radiation  component 
reaches  the  surface.  The  distribution  of  diffuse  radiation  over  the 
sky  dome  is  affected  by  variable  cloudiness  and  atmospheric 
clarity  conditions.  The  diffuse  solar  radiation  model  is  composed 
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Fig.  1.  Zenith  angle,  slope,  surface  azimuth  angle  and  solar  azimuth  angle  for  tilted  surface. 


Fig.  2.  Schematic  view  of  the  sky  showing  solar  radiance  distribution  component. 


of  three  parts:  isotropic  part  which  receives  solar  radiation 
uniformly  from  the  entire  sky  dome;  in  circumsolar  part,  diffusion 
takes  place  due  to  forward  scattering  of  solar  radiation;  horizon 
brightening  takes  place  in  clear  skies  and  near  the  horizon. 
Therefore,  the  angular  distribution  of  diffuse  radiation  is  a  function 
of  reflectance  (albedo)  of  ground.  A  high  albedo  value  results  in 
reflection  of  solar  radiation  back  to  the  sky  resulting  in  horizon 
brightening.  Thus,  the  solar  radiation  incident  on  tilted  surface  (/T) 
is  given  by  the  relation: 

h  =  ^T,b  T  ^T.d.iso  T  ^T,d,cs  "h  W,d,hz  T ^T.refl  (4) 

where  IT,b  is  the  beam  radiation,  IT,d,iso  is  the  isotropic  component, 
h,d,cs  is  the  circumsolar  component,  IT,d,hz  is  the  horizon  bright¬ 
ening  component  of  diffuse  radiation  and  ITrefi  is  reflected  radia¬ 
tion.  When  radiation  view  factor  is  taken  into  consideration,  (fT)  is 
as  follows: 

It  =  “ffd.iso^'c-s +fd,cs^£»  Tld./iz^c-liz  Tfpfc-g  (5) 

Pandey  and  Katiyar  [8]  investigated  different  diffuse  compo¬ 
nents  of  solar  radiation  at  different  tilt  angles  (15°,  30°,  45°  and 
60°)  at  Lucknow  (Latitude  26.75°,  Longitude  80.85°),  India.  The 
Circumsolar,  Isotropic  model,  Klucher  model  and  Hay's  anisotropic 
model  were  used  to  calculate  monthly  mean  daily  diffuse  radiation 
on  the  inclined  surface.  The  comparison  of  measured  and  esti¬ 
mated  monthly  mean  daily  diffuse  radiation  on  inclined  surface 
shows  that  the  Klucher  model  provides  a  good  estimation  of  the 
sky  diffuse  radiation  on  inclined  surface  for  northern  Indian  region 
conditions  with  mean  percentage  errors  of  less  than  7%.  David 
et  al.  [9]  used  four  (Hay,  Skartveit  and  Olseth,  Gueymard,  Perez) 


models  to  estimate  diffuse  radiation  on  tilted  surface  for  Reunion 
Island.  The  Perez  model  exhibits  best  performance  and  albedo 
constant  i.e.  0.2  provides  accurate  results.  Bortolini  et  al.  [10] 
developed  a  correlation  of  diffuse  fraction  in  terms  of  clearness 
index  to  estimate  diffuse  solar  radiation  for  European  sites.  The 
error  in  correlations  models  is  found  to  be  less  than  20%,  showing 
accurate  estimation  of  diffuse  radiation. 

A  number  of  models  have  been  proposed  by  different  research¬ 
ers  [11-25]  to  calculate  total  solar  radiation  on  tilted  surfaces  from 
global  radiation  on  horizontal  surface.  Temps  and  Coulson  [19] 
introduced  geometrical  terms  into  the  isotropic  model  to  take  into 
account  the  brightening  of  the  sky  in  the  region  of  the  sun  and  at 
the  horizon.  The  only  difference  among  these  models  is  in  the 
assessment  of  the  sky-diffuse  component.  Based  on  the  assump¬ 
tions  made,  the  estimation  models  can  be  classified  into  isotropic 
[12]  and  anisotropic  [13,15]  ones. 


2.4.1.  Isotropic  model 

This  model  approximates  the  condition  of  an  overcast  sky  and 
assumes  that  the  intensity  of  sky  diffuse  radiation  is  uniform  over 
the  sky  dome.  Liu  and  Jordan  [12]  modified  Hottel  and  Woertz's 
model  [11]  by  including  beam,  isotropic  diffuse  and  solar  radiation 
diffusely  reflected  from  the  ground  and  obtained  the  equation  for 
IT  as  follows: 


h  =  Ib^b+ld 


1-cos 

2 


(6) 


The  isotropic  sky  model  for  monthly  mean  solar  radiation  on 
unshaded  tilted  surface  is  considered  [18,19]  and  HT  is  given  as 
follows: 


Hr  =  HbRb+Hd 


1  +  cos  /? 


+  Hp 


(T-COSA) 


(7) 


For  the  surface  in  the  northern  hemisphere  sloped  toward  the 
equator,  i.e.  for  surfaces  with  z=0"  the  equation  for  Rb  is  given  as 
below. 


^  cos(</>-/f)cos  i5sinws  +  0r/180)wssin(</>-/?)sin  5 
6  —  cos  </>cos  i5sinws  +  (»/180)wssin  r/>sin  5 

where  ws  is  the  sunset  hour  angle  for  the  tilted  surface  for  the 
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mean  day  of  the  month  and  is  given  by 


cos_1(-tan</>tan5) 
cos-1  (— tan(0— /?)tan5) 


(9) 


For  surfaces  in  the  southern  hemisphere  sloped  toward  the 
equator,  with  >'=180°  Rb  is  given  as  follows. 


Rb 


cos(<^+/?)cos  (5sinWj+(a/180)wssin(^+/f)sin  <5 
cos  i/>cos  (5sinws+(ff/180)wssin  ij6sin  S 


(10) 


where  Fisa  modulating  function  given  by  F=  1  -(Hd/H)2.  When  the 
sky  is  overcast,  the  value  of  F  is  zero  and  relation  (19)  becomes 
isotropic.  When  value  of  F  tends  to  one,  Eq.  (19)  can  be  used  for  the 
clear  sky  conditions  as  given  by  Temps  and  Coulson  [19].  Hay’s 
anisotropic  model  [27]  consists  of  an  isotropic  as  well  as  circum¬ 
solar  component.  However,  horizon  brightening  as  incorporated 
by  Klucher  [15]  is  not  taken  into  consideration.  Therefore,  Hay 
presented  the  following  relation  for  the  daily  sky  diffuse  radiation 
on  tilted  surface: 


cos-1(-tan</>tan5) 
cos-1  (-tan(</>+/?)tan<5) 


(ID 


Hs  =  Hri 


^-^R6+l(l+cos/J)(l-^-^ 

'ext  ^  \  'ext 


(20) 


where  declination  angle  <5  =  (23.45ff/180)sin(2^(284+n/365)). 

The  isotropic  diffuse  model  presents  a  good  fit  to  solar  radia¬ 
tion  data  under  overcast  skies  and  tends  to  underestimate  IT  under 
clear  and  partly  cloudy  conditions.  Klein  and  Theilacker  (KT 
method)  [25]  assumed  both  diffuse  and  ground  reflected  radiation 
as  isotropic  to  calculate  average  radiation  on  sloped  surfaces.  This 
method  showed  improved  results  over  the  isotropic  method  when 
compared  with  hourly  calculation  of  radiation  data. 


g=  COS  UP-P) 
dc os  ip 


'  rrW,  „ 

SlnWs_T80COSW* 


b  (nWs 

2  V180' 


(12) 


where  wf  =  cos-1(-tan(0  +  /?)tan  <5),  the  positive  sign  is  used  in  the 
northern  hemisphere  for  south-facing  surfaces  and  the  negative 
sign  for  the  southern  hemisphere  for  north-facing  surfaces 
d  =  sinws-(^ws/180)cosws,  a  =  0.409  +  0.501 6sin(ws-60), 

b  =  0.6609-0.4767sin(ws-60). 

The  developed  equations  in  this  method  are  valid  for  any 
azimuth  angle  for  all  latitudes.  If  j^O”  or  180°  the  times  of  sunrise 
and  sunset  will  not  be  symmetrical  at  solar  noon;  Rb  is  formulated 
as  follows  and  it  is  valid  in  case  when  the  sun  rises  or  sets  on  the 
surface  twice  each  day: 


(13) 


max(0,G(wss,wsr))  if  wss>wsr 

max(0,(G(Wss,-wsr)  +  G(wss,wsr)))  if  wss  >  wsr 


(14) 


|  wsr  |  =  min 


ws,  cos 


_i  AB+C\/A*-B2  +  c2 


Vtf+B2 

wsr  =  —  I wsr |  if(/l  >  0  and  B  >  0)  or  (A>B)  |wsr|  otherwisej 


|wss|  =  min 


Ws,  COS 


_1  AB-C  s/A^-B2  +  C2 


V^+B2 


(15) 

(16) 
(17) 


Temps  and  Coulson  [19]  incorporated  horizon  brightening  on 
clear  days  by  using  correction  factor  (1  +  sin3(/?/2))  to  the  isotropic 
diffuse  model.  By  adding  modulation  factor  (f)  Reindl  et  al.  [14] 
modify  the  Hay  and  Davies  model  and  predicted  ldT  as  follows: 


Id.T  =  Id 


(1  -A)^1+2°S^  +  (l  +/sin3  (If) ^  +AjRb 


(21) 


where  /  =  s/Ib/I  and  A,  is  an  anisotropic  index  as  a  function  of 
atmosphere  transmittance  for  beam  radiation.  The  HDKR  (Hay, 
Davies,  Klucher,  Reindel)  model  considered  beam  and  ground 
reflected  terms  and  predicted  as  follows: 


(22) 


This  model  gives  total  radiation  on  tilted  surface  toward  the 
equator  and  gives  results  close  to  the  measured  values. 

Perez  et  al.’s  model  [20]  considered  all  the  components  of 
diffuse  radiation  and  (Idj)  is  as  follows. 


I(i  T  =  Id 


(1-F,)[l±f^)+F,!+F2Si„„' 


(23) 


where  a  =  max(0,cos  6)\  b  =  max(cos  85,cosdz)  accounts  for  the 
angle  of  incidence  of  the  cone  of  circumsolar  radiation;  brightness 
coefficients  are  given  by  F!  =  max[0,(fn +/12A+(®BZ/180)/13)]; 
f2=(f2i+/22A  +  (ffdz/180)/23);  clearness  e  =  ((ld+h  J/ld  +  5.535  x 
10_66)|)/  (1 +  5.535  x  10_6df);  brightness  A  =  m(/d/lon)  where  m 
and  Ion  are  the  air  mass  and  incidence  beam  radiation  respectively. 
Therefore,  by  adding  beam,  diffuse  from  isotropic,  circumsolar, 
horizon  and  ground  reflected  radiation  IT  is  given  by 

It  =  W+Ul-F,)^”5^  +/dFl  °  +IdF2 sin li+lP(^°^j 

(24) 

This  model  gives  slightly  higher  value  of  IT  in  comparison  to  the 
measured  values.  This  model  is  recommended  for  surfaces  with  y 
far  from  0°  in  the  northern  hemisphere  or  180°  in  the  southern 
hemisphere. 


wss  =  ^  |wss|  if(71  >  0  and  B  >  0)  or  (A>B) 

- 1  wss  |  otherwise]  (18) 

where  A  =  cos/3  +  tan^  cos^  sin/?,  B  =  coswscos^ + tan<5  sin/?  cos^, 
C  =  (sin/lsin^/cosi/i). 

This  method  is  useful  for  surfaces  with  surface  azimuth  angles 
more  than  15°  from  the  south. 


3.  Optimum  tilt  angle  determination 

The  optimum  tilt  angle  determination  of  a  solar  panel  is 
fundamental  to  its  efficient  operation  because  incorrect  position¬ 
ing  leads  to  loss  of  potential  solar  power.  The  optimum  tilt  angle 
calculations  are  based  on  maximizing  the  solar  radiation  falling  on 
a  sloped  surface  using  different  optimization  techniques. 


2.4.2.  Anisotropic  model 

Klucher  [15]  modified  the  Temps  and  Coulson  model  [19]  by 
incorporating  conditions  of  cloudy  skies  given  by  the  following 
relation  [26]: 

Hs  =  ^Hd(l  +cos/?)  ^1  +Fsin3(/?/2)J(l  +Fc  os2dsin3dz)  (19) 


3.3.  Total  solar  radiation  on  tilted  surface 

The  monthly  average  daily  total  solar  radiation  on  a  tilted 
surface  (Ffr)  is  the  sum  of  the  direct  or  beam  solar  radiation  (HB), 
diffuse  radiation  (HD),  and  ground  reflected  radiation  (HR)  on  the 
tilted  surface.  Thus,  for  a  surface  tilted  at  an  angle  from  the 
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horizontal,  the  incident  total  radiation  is  given  by  the  relation: 

Ht  =  Hb  +  Hd  +  Hr  (25) 

The  daily  beam  radiation  received  on  the  tilted  surface  may  be 
expressed  as 

HB  =  (Hg-Hd)Rb  (26) 

where  Hg  and  Hd  are  the  monthly  mean  daily  global  and  diffuse 
radiations  on  a  horizontal  surface.  The  Rb  for  the  surface  in  the 
northern  and  southern  hemispheres  sloped  toward  the  equator 
are  given  by  relations  (8)  and  (10). 

The  methods  to  estimate  the  ratio  of  diffuse  solar  radiation  on  a 
tilted  surface  to  that  of  a  horizontal  are  classified  as  the  isotropic 
and  anisotropic  models.  The  sky-diffuse  radiation  can  be 
expressed  as 

HD  =  RdHd,  (27) 

where  Rd  is  given  by  different  authors  as  shown  in  Table  1  [26]. 

As  a  consequence,  the  total  solar  radiation  on  a  tilted  surface 
can  thus  be  expressed  as 

HT  =  (Hg-Hd)Rb  +  Hgp1^°SP+HdRd  (28) 

Therefore,  total  solar  radiation  falling  on  tilted  surface  is 
computed  by  varying  tilt  angle  (/?)  from  0  to  90°.  The  optimum 
tilt  angle  is  taken  at  which  solar  radiation  on  the  tilted  surface  HT 
becomes  maximum. 


3.2.  Tilt  angle-latitude  relations 

A  number  of  authors  have  determined  the  optimum  angle- 
latitude  (<p)  relations,  which  are  generally  used  for  determining  the 
optimum  tilt  angles  used  as  thumb  rules  by  solar  energy  system 
installers  for  many  locations  (Table  2)  [42].  These  methods  give 
simple  and  approximate  values  for  tilt  angle;  however,  the 
accurate  optimum  angle  for  a  particular  site  should  be  determined 
using  the  measured  values  of  solar  radiation  for  the  site. 

The  latitude  based  optimum  tilt  angle  relations  do  not  take  into 
account  the  factors  that  affect  solar  radiation  like  altitude  above 
the  sea  level  and  cloud  coverage  patterns,  so  these  methods  do  not 
give  accurate  results  for  the  cloudy  regions  like  Northern  Europe. 
It  is  difficult  to  determine  their  accuracy  for  different  sites,  but 
these  relations  are  mostly  used  by  PV  industry  for  solar  system 
installation.  The  NREL's  PV  Watts  calculator— the  free  online  tool, 
[41  ],  can  also  be  used  with  site  specific  meteorological  data  to  find 
solar  energy  yield  at  any  tilt  and  the  optimum  tilt  is  taken  at  which 
solar  energy  yield  is  maximum.  Charles  Landau  recommended  that 
the  tilt  angles  in  winter  seasons  are  about  5°  steeper  than  the 
commonly  recommended  tilt  angle  [41].  In  winter  season,  the 
solar  radiation  is  mostly  incident  at  midday  so  panel  should  be 
mounted  at  solar  noon  to  collect  maximum  solar  energy  through¬ 
out  the  day.  In  summer  season,  the  sun  moves  lower  in  the  sky  in 
morning  and  evening  so  it  is  necessary  to  tilt  less  to  the  south  (or 
more  to  the  north)  to  collect  maximum  solar  radiation. 

The  calculation  of  tilt  angle  by  maximizing  solar  radiation  on 
sloped  surfaces  can  generally  be  done  for  sites  lying  between  5‘  N 


Table  1 

Values  of  parameter  Rd  for  different  sky  diffuse  models. 


Model 

Authors 

Mathematical  relations  for  Rd 

Isotropic 

Badescu  [28] 

Rd  ~  3+cf2W 

Isotropic 

Tian  et  al.  [29] 

Isotropic 

Koronakis  [30] 

Rd  =  J  [2+ COS/9] 

Isotropic 

Liu  and  Jordan  [24] 

Rt(=  1  +  COM 

Anisotropic 

Reindl  et  al.  [14] 

R*  =  +  (l-a)  (’T*)  (  (l  +  v/|sin3  (§)) 

Anisotropic 

Skarveit  and  Olseth  [17] 

Rd  —  0.511?/,  +  1 +2°s  ^  ^[sin/S  GtoJcos/?  *sin3(§)1 

where  Q  =  j  max  ^0,  ^0.3-2  J  j 

Anisotropic 

Steven  and  Unsworth  [31  ] 

Rd  =  0.51R/,+  1+“s^  i.2&.  [sm P  Gto )C0SP  'rS11’3©] 

Anisotropic 

Hay  [27] 

Ra^Rs  +  f+^X^) 

Anisotropic 

Klucher  [15] 

Rd  =  1(1+ cos  p)\\  -t-Fsin3 (>9/2)J  (1  +Fcos2<9sin30z) 

where  F  =  l-(/d//)2,cos  0  =  sin  <5sin(0-/?)+cos  <5cos(0-/?)cos  w. 

sin  0Z  =  ^-(sin  <5sin  0+cos  <5cos  0cos  w)2,w  =  cos-1  (-tan  0tan  S) 

Table  2 

Optimum  tilt  angle-latitude  angle  relations. 


References 

Recommended  tilt  angle 

Remarks 

Duffie  and  Beckman  [6] 

(0  +  15  )  + 15 

These  methods  are  suitable  for  calculating  approximate  tilt  angles 

Hey  wood  [32] 

0-10° 

for  different  latitudes.  The  minus  sign  is  used  for  summer 

Lunde  [33] 

0  + 15° 

season  and  the  positive  sign  for  winter  season. 

Chinnery  [34] 

0  +  10° 

Lof  and  Tybout  [35] 

0+(lO°-*3O°) 

Garg  [36] 

0+ 15",  (p— 15°,  0.90 

Hottel  [37] 

(p+ 20° 

Kern  and  Harris  [38] 

0+10° 

Yellot  [39] 

0+20° 

Elminir  et  al.  [40] 

(0  +  15  )  + 15 
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and  40  N  latitudes  with  sunny  climates  where  beam  radiation 
component  dominates.  The  areas  located  beyond  45°N,  like 
countries  in  Northern  Europe  namely  Ireland,  England,  Norway 
and  Sweden,  receive  maximum  amount  of  diffuse  radiation  due  to 
frequent  heavy  cloud  cover.  Therefore,  the  methodologies  of 
choosing  tilt  angle  by  maximizing  the  beam  radiation  are  unsui¬ 
table,  as  other  atmospheric  conditions  are  not  considered. 

Armstrong  and  Hurley  [42]  developed  a  methodology  to 
determine  optimum  tilt  angle  for  locations  with  frequently  over¬ 
cast  skies  using  monthly  sunshine  duration  data  and  hourly  cloud 
observations.  Under  cloudy  skies,  it  is  important  to  differentiate 
between  direct  and  diffuse  radiation  for  a  particular  site  to 
calculate  optimum  tilt  angle  so  the  Perez  model  is  useful  to 
calculate  diffuse  radiation  falling  on  the  solar  panel.  The  tilt  angle 
is  changed  from  0°  to  90°  in  steps  of  1°  and  is  selected  as  the  one 
that  maximizes  the  incident  solar  radiation  on  solar  panel  by 
taking  into  account  the  frequencies  of  cloudy  skies.  The  optimum 
tilt  angle  of  grid  connected  and  standalone  PV  system  is  selected 
that  matches  the  available  solar  radiation  with  the  load  demand. 

3.3.  Tilt  angle  optimization  using  CA,  SA  and  PSO  techniques 

Optimization  is  a  method  of  finding  the  conditions  that  give 
maximum  or  minimum  value  of  the  objective  function  to  obtain 
the  accurate  results  under  specified  conditions.  It  is  useful  in  the 
design,  construction  and  maintenance  of  any  engineering  system 
[43].  For  finding  optimum  tilt  angle,  solar  radiation  on  the  tilted 
surface  (Eq.  (28))  is  taken  as  an  objective  function  which  is  solved 
using  different  optimization  techniques  like  Genetic  Algorithm 
(GA),  Simulated  Annealing  (SA)  technique  and  particle  swarm 
optimization  (PSO). 

Genetic  algorithm  (GA)  is  suitable  for  optimization  problems 
which  include  complex  nonlinear  variables,  as  it  can  find  the 
global  optimum  solution  with  a  high  probability  [44,45].  A 
population  of  points  is  used  for  starting  the  GA  instead  of  a  single 
design  point  [46].  It  involves  principles  of  natural  genetics  and 
natural  selection.  The  natural  genetics  are  reproduction,  crossover 
and  mutation  which  are  used  in  the  genetic  search  procedure  [47]. 
Talebizadeha  et  al.  [48]  used  GA  to  calculate  hourly,  daily,  monthly, 
seasonally  and  yearly  optimum  tilt  angle  for  Iran  and  showed  that 
the  optimum  hourly  surface  azimuth  angle  is  not  zero  and 
optimum  tilt  angles  of  photovoltaic  panels  and  solar  collector  are 
found  to  be  the  same.  The  solar  energy  gain  at  daily,  monthly 
optimum  tilt  angle  is  found  to  be  the  same  but  energy  gain  is  more 
at  hourly  tilt  angle.  Therefore,  the  solar  tracker  is  beneficial  for 
hourly  variation  of  tilt  angle.  Congradac  et  al.  [49]  used  GA  and 
fuzzy  logic  process  to  find  out  optimum  blind  tilt  angle  (angle 
rotated  in  anticlockwise  and  clockwise  direction)  for  maintaining 
accurate  brightness  of  the  room.  This  process  is  useful  in  main¬ 
taining  user's  comfort  and  saving  energy. 

Simulated  Annealing  (SA)  derives  its  name  from  the  simulation  of 
thermal  annealing  of  critically  heated  solids  and  is  used  to  find  the 
global  optimum  with  a  high  probability  of  objective  functions  which 
contain  numerous  local  minima.  Chen  et  al.  [50]  used  SA  for 
calculating  the  optimum  installation  angle  for  fixed  solar-cell  panels. 

Particle  Swarm  Optimization  (PSO)  is  a  stochastic  technique  for 
exploring  the  search  space  for  optimization  [51].  It  is  achieved  by 
particles  in  multidimensional  space  that  have  a  position  and  a 
velocity  [52].  In  PSO  particles  are  flying  through  hyperspace  and 
adjust  their  own  best  position  i.e.  the  smallest  objective  function 
values.  Chang  [53]  used  the  varying  inertia  weight  methods  [54— 
58]  and  proposed  a  particle-swarm  optimization  method  with 
nonlinear  time-varying  evolution  (PSO— NTVE)  to  determine  the 
tilt  angle  of  PV  modules  for  maximum  output  electrical  energy  of 
the  modules  in  Taiwan.  The  yearly  optimal  angle  for  Taipei  area  is 
18.16°  and  17.3°,  16.15°,  15.79°,  15.17°,  17.16°,  15.94  for  Taichung, 


Tainan,  Kaosiung,  Hengchung,  Hualian,  Taitung  respectively.  The 
PSO— NTVE  is  faster  than  the  other  three  PSO  methods  and  the  GA 
method  for  achieving  an  optimum  solution.  Therefore,  these 
optimization  techniques  will  give  better  results  than  the  other 
methods  (latitude  based  method,  maximizing  radiation  by  chan¬ 
ging  /I  from  0°  to  90"  at  fixed  steps). 

3.4.  Tilt  angle  optimization  using  ANN  techniques 

Artificial  Neural  Network  techniques  [ANN]  are  excellent  tools  for 
research  as  these  are  able  to  solve  non-linear  function  approximation, 
data  classification,  clustering  and  simulation  [59],  These  are  used  in 
different  fields  of  science  and  technology  [60-63]  and  also  for  the  solar 
radiation  predictions  [64].  Chang  [65]  used  sequential  neural-network 
approximation  and  orthogonal  arrays  (SNAOA)  to  determine  the 
optimum  tilt  angle  for  maximum  output  power  energy  of  PV  modules 
in  seven  areas  of  Taiwan.  The  back  propagation  neural  network  is  used 
to  search  optimum  tilt  angle  design.  The  maximum  and  minimum 
values  of  optimum  tilt  angle  and  voltage  are  used  as  constraint  to 
maximize  power  on  PV  module.  The  annual  optimal  angle  is  found  to 
be  23.25°  for  Taipei  site.  Mehleri  et  al.  [66]  used  Radial  Basis  Function 
neural  network  (RBFNN)  to  predict  solar  radiation  on  tilted  surface 
using  tilt  angle  and  orientation  as  input  parameters.  Fuzzy  algorithm  is 
used  for  training  RBFNN  and  it  is  found  that  changing  tilt  angle 
seasonally  improves  performance  and  uniformity  of  PV  array  output 
power.  Mehleri  et  al.  [67]  incorporates  global  solar  irradiance  on  a 
horizontal  surface,  extraterrestrial  radiation,  solar  zenith  angle  and 
solar  incidence  angle  on  a  tilted  plane  as  inputs  to  RBFNN  for 
estimating  global  solar  irradiance  on  inclined  surfaces  in  Athens.  The 
absolute  fraction  of  variance  (R2)  is  96%,  proving  the  estimation  by  the 
model  accurate.  Notton  et  al.  [68]  developed  three  ANN  models  for 
Mediterranean  site  of  Ajaccio,  France,  to  estimate  hourly  global 
irradiation  on  the  tilted  plane.  The  first  and  second  ANN  models 
estimate  hourly  global  irradiation  on  45°,  60  tilted  planes  respectively 
whereas  the  third  ANN  model  estimates  both  on  45°,  60°  tilted  planes. 
The  ANN  models  incorporate  declination  angle,  hour,  zenith  angle, 
hourly  extraterrestrial  horizontal  irradiation,  hourly  global  irradiation 
as  inputs  and  hourly  global  irradiation  on  45°,  60  tilted  planes  as 
output  parameters.  The  (R2)  in  the  first,  second,  third  ANN  models  are 
99.79%,  99.82%  and  99.70%  respectively.  Celik  and  Muneer  [69]  used 
generalized  regression  neural  networks  (GRNN)  to  predict  solar 
radiation  on  the  tilted  surface  for  Iskenderun,  Turkey.  The  GRNN 
utilize  input  parameters  as  global  solar  irradiation  on  horizontal 
surface,  declination  and  hour  angles.  The  R2  and  mean  absolue 
percentage  error  are  98.7%  and  14.9  Wh/m2  respectively.  Chatterjee 
and  Keyhani  [70]  used  14  inputs  (latitude,  ground  reflectivity  and  12 
months  irradiance  value)  to  estimate  optimum  tilt  and  total  irradiance 
on  tilted  surface  by  ANN.  The  activation  function  in  hidden,  output 
layers  is  hyperbolic  tangents,  and  linear  and  Levenberg-Marquardt 
algorithm  is  used  for  training.  The  difference  between  predicted  and 
analytically  determined  optimum  tilt  angle  is  found  to  be  3°  showing 
accurate  estimation  by  the  ANN  techniques  Therefore,  the  ANN 
techniques  reduce  the  computation  part  of  optimum  tilt  angle 
calculations,  and  estimate  optimum  angle  and  tilted  surface  solar 
radiation  at  different  sites  with  better  accuracy. 

4.  Overview  of  optimum  tilt  angle  determination  methods 

Many  authors  have  determined  optimum  tilt  angles  analytically 
and  experimentally  by  searching  for  the  maximum  total  solar 
radiation  falling  on  the  PV  modules  for  a  number  of  locations. 
Some  of  the  relevant  highlights  are  given  in  this  section.  Asl- 
Soleimani  et  al.  [71]  determined  optimum  tilt  angle  for  Tehran 
using  an  experimental  setup  of  five  11 -Wp  PV  mono-crystalline 
modules  by  measuring  the  maximum  energy  output  for  a  period  of 
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Table  3 

Summary  of  optimum  tilt  angle  correlations. 


Authors  Description 


The  correlation  of  optimum  tilt  angle  popt  in  terms  of  clearness  index  ( KT ),  latitude  (0),  day  number  (D)  is  given 
by  'Popt  =  (6-4.8Kr  +  0 .86Kj27  0  +  0 . 002 1 02 ) 

+(31K?'37  +  0.094Jf?46^+0.00063KJ°17^2)cos^||5(D+11.5)j 

The  optimum  tilt  angle  (popt)  as  a  function  of  clearness  index  ( KT ),  declination  angle  (<5)  for  Gebze  in  Turkey  are  as  follows,  0opt=  33.24-1.315; 
popt  =  22. 09  +  25. 79/Cj-l. 495;  /9opt=35.15-1.37t5-0.007<52;  /?opt=35.15-1.39<5-0.007<52-4.26  x  10-5<53.  The  variation  of  root  mean  square  error 
(RMSE)  is  from  1.09  to  1.74  and  R2  from  99.40  to  99.77. 

The  monthly  popt  in  terms  of  latitude  (0),  Julian  Day  (M)  is  given  by  (i)  From  January  to  March  (1<M<3) 

popt=  60.0001 2  +  1.49986M + 3.49996M2 + (0-30)  x  (0.7901 +0.01794M+0.0165M2).  (ii)  From  April  to  June  (3<M<6)  popt= 216.0786 
-72.O3219M+6.OO312M2  +  (0-4O)  x  (1.07515  +  0.11244M-0.03749M2).  (iii)  From  July  to  September  (6<M<9)  popt= 29.118311 
-2O.52981M+2.5O186M2  +  (0-5O)  x  (-11.17256+2.70569M-0.15035M2).  (iv)  From  October  to  December  (9<M<12) 

popt=  -441.2385  +  84.54332M— 3.50196M2  +  (0—40)  x  (4.2137-0.54834M+0.0223M2).  Providing  that  if  popt  <  0  then  popt=0.  The  design  of  solar 
collector  at  negative  tilt  angle  is  significant  between  latitude  0°  and  +  25°. 

The  monthly  optimum  tilt  angles  of  the  south  facing  collectors  as  a  sinusoidal  function  of  latitude  (0),  and  day  of  the  year  D  over  Turkey  is  given 
by/?opf  =  25. 521438+26. 838291cos(-0. 0178440+1. 013901D+7. 527742).  The  R2,  RMSE  of  correlation  equation  are  98.8%  and  2.06%  respectively. 

The  daily  (popt,d),  monthly  ( popt,m )  optimum  tilt  angle  for  Syria  are  as  follows:  popUd  =  0-tan_1  [s  J“ss)  tan(g)] 

i  (  En  -  n?  ((24/jr)/sc[l  +0.034cos(2;rn/365)sin(<5)Wss)  \  _ 

pnntm  =cb- tan-1  n  =  \ — - — — - : -  where  n i,  n2  are  the  first,  last  days  of  the  mth  month  as  counted  from  January.  The 

i  opt.m  r  \En  =  n]  ((24/^)/sc[l  +0.034cos(2^n/365)cos(^)sin(wSs))  J  "  z  J  J  J 

yearly  optimum  tilt  angle  is  30.56°  and  energy  gains  are  28%,  26%  at  monthly,  season  tilt  angle. 

Moghadam  et  al.  [88]  The  annual  optimum  p0pt=  0.197 +0.321.  The  tilt  angle  is  negative  for  days  (n  =  50-130)  in  Zahedan  city  which  means  the  face  of  collector  is  in 
north  direction  during  these  days.  The  two  time  (20  Mar  to  21  Sep,  22  Sep  to  20  Mar)  adjustment  of  tilt  angle  increase  the  energy  gain  by  8%. 

Agha  and  Sbita  [89]  The  optimum  tilt  angle  is  selected  that  gives  best  matching  load  curve  with  solar  irradiation  curve  and  is  given  by  following  optimization  factor 
{EoptiP)),  in  terms  of  monthly  average  daily  total  solar  irradiation  (Havg/ip)),  curve  standard  deviation  ( SDEV(p )).  Fopt{p)=HaV^p)ISDEV(p) 

Talebizadeh  et  al.  [90]  The  seasonal  optimum  slope  angle  developed  by  correlation  is  popt=  1.0730+ 10.3  in  winter  season,  popt= 0.48850-10.27  in  spring  season, 

Popt= 0.26370+4.961  in  summer  season,  popt= 0.89660 +23.81  in  autumn  season  and  popt= 0.68040 +7.203  is  yearly  optimum  tilt  angle.  The 
optimum  daily,  monthly,  seasonally  and  yearly  azimuth  angle  is  zero  for  Iran.  The  energy  gain  is  maximum  at  daily  optimum  tilt  angle  and  solar 
tracking  system  is  not  economical  due  to  high  cost. 


Elsayed  [83] 

Tiris  and  Tiris  [84] 

Elkassaby  [85] 

Ertekin  et  al.  [86] 
Skeiker  [87] 


one  year  and  emphasized  the  need  for  validating  the  analytically 
determined  optimal  tilt  angles  experimentally  also.  Li  et  al.  [72] 
proposed  a  sunshine  hour  data  based  model  to  calculate  solar 
radiation  on  the  tilted  surface  with  different  orientations.  The 
maximum  error  in  this  model  is  found  to  be  less  than  5.2%  so  it  can 
be  used  to  determine  the  optimum  tilt  angle  for  sites  where  only 
sunshine  hour  data  are  measured.  Chang  [73]  used  orthogonal 
array  experiment  technique  and  an  ant  direction  hybrid  differen¬ 
tial  evolution  algorithm  (ADHDEOA)  for  tilt  angle  of  PV  modules  in 
Taipei  area.  An  ant  colony  method  is  used  to  find  the  optimum  tilt 
angle  at  which  power  is  maximum.  The  annual  optimal  angle  for 
Taipei  area  is  determined  (Taichung  17.3°,  Tainan  16.15°,  Kaosiung 
15.79°,  Hengchung  15.17°,  Hualian  17.16°,  Taitung  15.94°).  The 
computer  simulated  optimum  angles  are  also  verified  experimen¬ 
tally.  The  ADHDEOA  is  found  to  be  better  than  GA  method  in 
computation  and  is  much  effective  in  optimum  tilt  angle 
determination. 

Beringer  et  al.  [74]  carried  out  the  experimental  measurements 
for  monthly  optimum  tilt  angles  for  cloudy  mid-latitude  location 
Hannover,  Germany,  using  eight  multi-crystalline  silicium  solar 
photovoltaic  panels  with  south  orientation  mounted  at  different 
tilt  angles  from  0°  to  70°  in  steps  of  10°.  The  maximum  power  is 
found  to  be  at  50-70  tilt  angles  in  winter  and  0-30°  in  summer 
months.  The  yearly  performance  difference  is  found  to  be  less  than 
6%  showing  minor  differences  at  various  tilt  angles.  However, 
further  follow-up  studies  to  determine  optimal  tilt  angles  for  such 
locations  with  cloudy  conditions  can  further  be  taken  up  using  the 
optimization  techniques.  Pour  et  al.  [75]  used  isotropic  Liu  and 
Jordan  model  [12]  to  determine  the  optimum  tilt  angle  and  the 
non-isotropic  Klein  [23]  and  Hay  [13]  models,  which  consider  the 
effects  of  the  azimuth  angle,  are  used  to  determine  the  optimum 
tilt  and  azimuth  angles  for  Isfahan,  Iran  [lat  32  N,  Long  51°  E[.  The 
optimum  tilt  angles  determined  using  the  Liu  and  Jordan  model 


for  south  orientation  shows  that  the  yearly  optimum  tilt  angle  is 
nearly  equal  to  the  latitude  of  the  site,  whereas  the  results  using 
the  non-isotropic  models  i.e.  the  Klein  and  Hay  models  show  that 
the  panels  at  optimum  tilt  and  azimuth  angles  receive  higher  solar 
energy  compared  with  the  fixed  solar  panel  which  indicates  that 
the  determination  of  both  optimum  tilt  and  azimuth  angles  is 
important  for  maximizing  incident  solar  radiation  at  a  location. 

A  number  of  studies  have  been  carried  out  to  maximize  power 
output  and  cooling  load  reduction  of  building  integrated  photo¬ 
voltaic  [BIPV]  as  such  the  optimum  tilt  angle  determination  for 
BIPV  design  and  development  is  another  important  area  of  PV 
applications.  Elhassan  et  al.  [76]  determined  optimum  tilt  angles 
for  building  integrated  photovoltaic  design  and  applications  at 
Kuala  Lumpur,  Malaysia  [Lat  2°  30'  N  Long  112°  30'  E  ],  by  using  a 
four  PV  module  experimental  setup  inclined  in  the  North,  South, 
East  and  West  directions  to  determine  the  optimum  tilt  angles  for 
these  directions  as  solar  radiation  received  during  the  day  varies 
for  each  direction.  The  optimum  tilt  angle  for  this  location  was 
found  to  be  nearly  equal  to  the  latitude  of  the  location.  Sunderan 
et  al.  [77]  investigated  that  the  optimum  tilt  angle  and  orientation 
enhance  the  power  generation  of  Standalone  Photovoltaic  Elec¬ 
tricity  Generation  Systems  (SPVEGS)  in  Ipoh,  Malaysia.  The  orien¬ 
tation  of  PV  module  is  North  from  April  to  August  and  south  for 
the  other  months.  The  energy  gains  of  6.4%  and  6.1%  are  achieved 
for  monthly  and  yearly  optimum  tilt  angles  respectively.  Lucio 
et  al.  [78]  have  given  an  algorithm  to  evaluate  optimum  tilt  angle 
for  obtaining  minimum  loss  of  load  probability  and  optimum 
design  of  stand-alone  photovoltaic  systems  in  Europe.  Asowata 
et  al.  [79]  used  a  photovoltaic  experimental  setup  at  Vaal  Triangle, 
South  Africa,  lying  in  the  southern  hemisphere  [  Lat.  29°  00'  S  and 
Long  24"  00'  E  ]  to  validate  the  optimum  tilt  angles  16°,  26°  and  36 
°  determined  analytically  for  winters  for  this  location.  Sun  et  al. 
[80]  determined  the  optimum  tilt  angles  for  the  shading-type  BIPV 
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Table  4 

Optimum  tilt  angles  for  solar  panels  for  different  locations. 


Study 

Location 

Monthly  optimum 
tilt  angle 

Yearly  optimum  tilt  angle 

Remarks 

Tsalides  and  Thanailakis  [98] 

Greece 

±  60°  (Greater  than  latitude) 

-  Optimum  tilt  angle  decreases  with  large 
azimuth  angles. 

-  Solar  radiation  uniformity  is  considered  for 
yearly  optimum  tilt  angle  calculation  resulting 
in  continuous  supply  of  power  to  load. 

Benghanem  [99] 

Madinah,  Saudi 
Arabia 

Winter  37°  Summer 
12° 

23.5°  equal  to  latitude 

-  At  yearly  optimum  tilt  angle,  the  energy  loss  is 
8%  in  comparison  to  monthly  optimum 

tilt  angle. 

-  For  more  accurate  tilt  angle  anisotropic  models, 
optimization  techniques  need  to  be  considered. 

Rowlands  et  al.  [100] 

Ontario  Canada 

Equal  to  latitude 

-  Azimuth  angle  is  close  to  south. 

-  PV  module  power  in  outdoor  conditions  should 
be  investigated. 

Kaldellis  and  Zafirakis  [101] 

Athens 

15°  with  (  ±  25°) 
deviation  for 

summer 

-  Performance  difference  of  PV  at  fixed  and 
variable  optimum  tilt  angle  is  27.3%. 

Calabroa  [102] 

USA  and  Europe 

pop[=(f)-( 26°,  27°,  28°)  where  $  varies  from  36° 
to  46 

-  Semi-fixed  solar  panels  are  beneficial  due  to 
different  seasonal  angles. 

Agarwal  et  al.  [103] 

Delhi  and  Nandha 

Tilt  angle  is 

Calculated  angle  by  Liu  and  Jorden,  Reindel, 

-  Incident  solar  radiation  is  maximum  at  monthly 

in  Haryana 

minimum  in 
summer  and 
maximum  in  winter 

season. 

Hay  and  Badescu  models  is  28.29°,  29.41  , 

29.08 , 27.58  for  Delhi  and  30.61°,  31.66°, 
31.20°,  29.58°  for  Nandha. 

tilt  angles  and  minimum  at  yearly  optimum 
tilt  angle. 

Zhao  et  al.  [104] 

Singapore 

Daily,  monthly 
average  varies  from 
0.1°  to  25.9°,  0.1°  to 
28.7°. 

-  Energy  (kWh/m2-year)  received  at  0°,  monthly, 
daily  tilt  are  1264.78, 1297.32, 1298.86  proving 
monthly  tilt  angle  change,  is  beneficial. 

Daut  et  al.  [105] 

Perlis,  Northern 
Malaysia 

-17.160°  to  29.740° 

(  ±  means  PV  facing 
is  north,  south) 

-  Popt=(i,~3- 

-  Changing  PV  orientation  and  tilt  is  necessary. 

-  Monthly  beam  solar  irradiance  (968.36  W/m2) 
can  be  increased  at  optimum  tilt  angle  so  good 
scope  of  PV  generation. 

Khatib  et  al.  [106] 

Kuala  Lumpur. 

Johor  Bharu,  Ipoh, 
Kuching ,  Alor  Setar 
in  Malaysia 

Maximum  in 
summer  and 
minimum  in  winter 

-  Monthly  optimum  tilt  angle  is  recommended. 

Kacira  et  al.  [107] 

Sanliurfa,  Turkey 

13  and  61  for  June 
and  December 
months 
respectively. 

-  Solar  radiation  gain  at  monthly  optimum  tilts 
angles  is  more  than  seasonal  and  yearly  tilt 
angles.  In  two-axis  solar  tracking  system,  there 
is  a  gain  of  29.3%  in  solar  radiation  resulting  in 
34.6%  gain  in  generated  power. 

-  Optimization  techniques  can  be  used  for  better 
accuracy. 

Yakup  and  Malik  [108] 

Brunei  Darussalam 

Sep  1.6°,  Dec  32.3° 

33° 

-  Monthly  tilt  angle  is  recommended. 

-  Correlation  of  solar  noon  optimum  tilt  angle  in 
terms  of  latitude  are  given,  which  are  useful  to 
collect  maximum  radiation  in  winter 

Ahmad  and  Tiwari  [109] 

New  Delhi 

Winter  47.5° 
(0+90°),  summer 
13°  (0-60°) 

Equal  to  the  latitude 

-  The  loss  of  energy  is  1  %  and  15  %  when  tilt 
angle  is  adjusted  seasonal  and  yearly  instead  of 
monthly  optimum  tilt  angle. 

-  Using  measured  diffuse  radiation  in  Anisotropic 
models  can  improve  the  results 

Gunerhan  and  Hepbasli  [110] 

Izmir,  Turkey 

Summer 

Popt=<P+ 15°  winter 
Popt  =  (P~  15 

Equal  to  latitude 

-  Extraterrestrial  radiation  is  used  for  tit  angle 
calculation  resulting  in  inaccuracy. 

Ulgen  [111] 

Izmir,  Turkey 

June  0°,  December 

61 

30.30° 

-  Total  solar  radiation  at  0°,  61°,  30.30°  tilt  angles 
are  27.07  MJ/m2/day,  10.81  MJ/m2/day  and 
6397.78  MJ/m2/year. 

Elhab  et  al.  [112] 

Kuala  Lumpur, 
Malaysia. 

p=  10° 

-  Cooper's  equation,  Visual  Basic  Application 
programming  are  used 

Xianping  [113] 

Changsha 

19.22°  i.e.  10°  less  than  latitude. 

-  Anisotropic  Hay's  model,  p  of  0.14  (Magallanite) 
is  used  to  calculate  optimum  tilt  angle. 
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claddings  at  different  orientations  for  Hong  Kong  [lat  22.25  N  Long 
114.1667°S].  The  maximum  electricity  generation  per  unit  PV  area 
is  found  when  the  PV  modules  are  installed  on  south  facades  at 
the  tilt  angle  of  10°,  thus  increasing  the  energy  gain  significantly. 

Bojic  et  al.  [81]  determined  optimum  tilt  angles  for  PV  systems 
that  are  located  in  four  towns  (Les  Avirons,  Petite-France,  Saint- 
Benoit.and  Piton  Saint-Leu)  in  Reunion  Island,  France.  The  Hooke- 
Jeeves  algorithm  is  used  to  find  optimum  tilt  angle  at  which 
electrical  energy  generation  becomes  maximum  and  found  a 
discrepancy  between  latitude  and  optimum  tilt  angles,  showing 
the  need  for  determining  optimum  tilt  angles  for  every  location. 
Siraki  and  Pillay  [82]  proposed  a  modified  HDKR  (Hay,  Davies, 
Klucher,  Reindl  model)  anisotropic  sky  model  which  also  considers 
the  effects  of  the  adjacent  buildings  in  urban  locations  to  calculate 
the  optimum  angle.  The  optimum  angles  for  the  five  different 
latitudes  were  determined  and  the  tilt  angle  dependence  on 
latitude,  weather  conditions,  surrounding  obstacles  and  optimum 
azimuth  was  found.  The  optimum  tilt  angle  is  found  to  be  close  to 
the  location  latitude  for  small  latitude,  while  for  higher  latitude 
the  optimum  angle  is  found  to  be  smaller  than  the  location 
latitude.  These  studies  indicate  that  for  maximum  solar  radiation 
capturing  both  orientation  and  optimum  tilt  angle  of  PV  array 
needs  to  be  determined  accurately  for  any  location. 

Several  authors  have  given  correlations  of  optimum  tilt  angle 
[Pope)  in  terms  of  different  parameters  for  the  calculation  of 
optimum  tilt  angles,  as  shown  in  Table  3. 

Hartely  et  al.  [91]  compared  Temps-Coulson's,  Klucher's  and 
Hay's  models  for  the  estimation  of  total  radiation  on  vertical 
planes  facing  north,  south,  east  and  west  to  find  the  best  model 
that  can  be  used  to  calculate  solar  irradiance  on  tilted  planes  in 
Valencia  Spain  [Lat  39.5  N,  Long  0.67°  W  ].  Hay's  model  was  found 
to  be  accurate  in  determining  hourly  variations  of  optimum  tilt 
angles  for  south  facing  collector  but  the  average  monthly  optimum 
tilt  angle  was  not  given  accurately.  Tang  and  Wu  [92]  developed  a 
method  to  estimate  the  optimal  tilt  angles  of  a  solar  collector 
based  on  the  monthly  horizontal  radiation  and  prepared  a  map 
showing  optimal  tilt  angles  for  the  year  for  south  facing  collectors 
in  China  and  used  Collares-Perara  and  Rabl  [93]  correlation  to 
estimate  diffuse  solar  radiation  and  Klein's  method  [23]  for  the 
calculation  of  Rb.  It  is  found  that  Collares-Perara  and  Rabl 
correlations  give  accurate  optimum  tilt  angle  except  for  low 
clearness  index  sites. 

Fahl  and  Ganapathisubbu  [94]  computed  optimum  tilt  angles 
for  Bangalore  [Lat  12.97"  N,  Long  77.56°  E]  under  various  tracking 
conditions.  The  optimum  tilt  angles  are  estimated  between  15° 
and  17°  for  south  facing  collectors.  The  incident  solar  radiation  is 
found  to  be  the  same  within  variation  of  +3°  in  optimum  tilt 
angle.  The  output  of  a  monthly  adjusted  collector  at  optimum 
angle  is  found  to  be  increased  by  10%  in  comparison  to  that 
received  at  horizontal  orientation.  But  the  output  is  found  to  be 
35%  more  for  collectors  with  tracking  system. 

Chang  [95]  calculated  the  optimal  tilt  angles  for  six  locations  in 
Taiwan  using  the  extraterrestrial,  predicted  global  radiation  model 
and  ten-year  measured  solar  radiation  data.  The  results  indicate 
that  the  tilt  angles  calculated  from  extraterrestrial,  predicted  solar 
radiation  are  latitude  dependent  whereas  those  calculated  using 
measured  data  vary  from  site  to  site  and  shows  that  the  optimal 
tilt  angle  is  accurately  determined  using  measured  data  for  a 
location  with  cloudy  or  pollutant  environment.  The  method  of 
maximizing  extraterrestrial  radiation  for  tilt  angle  calculation  does 
not  consider  attenuation  of  the  solar  radiation  so  it  will  not  give 
accurate  tilt  angles.  Shariah  et  al.  [96]  used  the  annual  solar 
fraction  as  an  indicator  to  find  the  optimum  tilt  angles  for  a 
thermosyphoning  solar  water  heater  in  the  northern  and  southern 
parts  of  Jordan.  The  optimum  tilt  angle  varies  from  4>  to  tp  ±  20°  for 
high  solar  fraction  system.  Tang  et  al.  [97]  developed  a 


mathematical  procedure  for  determining  optimum  tilt  angles 
based  on  maximum  annual  collectible  radiation  for  evacuated 
glass  tube  collectors  in  China  and  found  that  the  optimum  tilt 
angles  are  found  to  be  about  10°  less  than  the  latitude  for  sites 
with  latitudes  greater  than  30°. 

The  optimum  tilt  angle  study  highlights  for  different  locations 
by  the  other  authors  are  summarized  in  Table  4. 


5.  Conclusions 

In  this  paper,  optimum  solar  panel  tilt  angle  methods  using 
different  optimization  techniques  are  reviewed.  Based  on  the 
study,  the  main  conclusions  are  as  follows: 

•  Optimum  tilt  angle  determination  maximizes  the  incident  solar 
radiation,  resulting  in  optimum  sizing  of  solar  systems. 

•  The  optimum  tilt  angles  are  determined  by  changing  tilt  angle 
from  0°  to  90°  in  different  steps  and  using  GA,  SA,  PSO  and  ANN 
optimization  techniques. 

•  The  latitude  based  thumb  rules  mostly  used  by  PV  and  solar 
thermal  industry  may  be  suitable  for  certain  locations  but 
result  in  escalating  the  system  costs  and  oversizing  of  systems 
if  used  without  proper  analysis.  This  effect  will  be  more 
prominent  in  case  of  solar  power  plants  where  a  large  number 
of  solar  panels  are  used. 

•  The  energy  gain  is  found  to  be  more  at  monthly  optimal  tilt 
angle  than  seasonal  and  yearly  tilt  angles  for  any  location. 

•  The  ANN  models  can  further  be  investigated  using  different 
combinations  of  input  variables  such  as  horizontal  solar  radia¬ 
tion,  latitude,  extraterrestrial  radiation,  the  isotropic  and  aniso¬ 
tropic  models,  diffuse  and  beam  radiation  for  estimating 
optimum  tilt  angle. 

•  For  better  accuracy  of  optimum  tilt  angle  calculations  different 
anisotropic  models  and  optimization  techniques  can  be  tested 
at  different  sites. 

•  For  urban  areas,  the  obstacles  affecting  the  solar  radiation 
should  also  be  considered  for  computing  optimum  tilt  angles. 

The  results  of  different  optimization  and  ANN  techniques  can 
be  compared  to  identify  the  best  methods  for  further  research. 
Although,  a  number  of  studies  have  been  carried  out  using 
predicted  solar  radiation  data,  extraterrestrial  radiation  data  but 
the  determination  of  these  angles  using  measured  global,  diffuse 
and  beam  solar  radiation  values  will  result  in  better  accuracy. 
However,  this  may  require  elaborate  reliable  solar  radiation 
measuring  equipments  at  the  location.  It  will  be  interesting  to 
study  the  effect  of  tilt  angles  on  annual  energy  production  and 
economics  of  solar  power  plants  installed  at  various  locations 
worldwide.  This  will  provide  valuable  inputs  to  researcher  and 
solar  industry  to  install  efficient  and  cost  effective  solar  systems  at 
different  sites  using  optimal  tilt  angles  for  each  location 
worldwide. 
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